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CONVEX INTEGRATION AND INFINITELY MANY WEAK 
SOLUTIONS TO THE PERONA-MALIK EQUATION 
IN ALL DIMENSIONS 


SEONGHAK KIM AND BAISHENG YAN 


Abstract. We prove that for all smooth nonconstant initial data the initial- 
Neumann boundary value problem for the Perona-Malik equation in image process¬ 
ing possesses infinitely many Lipschitz weak solutions on smooth bounded convex 
domains in all dimensions. Such existence results have not been known except for 
the one-dimensional problems. Our approach is motivated by reformulating the 
Perona-Malik equation as a nonhomogeneous partial differential inclusion with lin¬ 
ear constraint and uncontrollable components of gradient. We establish a general 
existence result by a suitable Baire’s category method under a pivotal density hy¬ 
pothesis. We finally fulfill this density hypothesis by convex integration based on 
certain approximations from an explicit formula of lamination convex hull of some 
matrix set involved. 


1. Introduction 


In this paper, we study the initial and Neumann boundary value problem: 



( 1 . 1 ) 


where hi C R n is a smooth bounded convex domain, T > 0 is a given number, 
u = u(x,t ) is the unknown function with u t denoting its time-derivative and Du = 
(u xi , • • ■ , u Xn ) its spatial gradient, n is outer unit normal on dD, and u 0 (x) is a given 
smooth function satisfying 


( 1 . 2 ) 


Duq ^ 0 in H, duo/dn = 0 on dfl. 


Problem (11.11) . especially when n = 2, is a famous Perona-Malik model in image 
processing introduced by Perona and Malik [27] for denoising and edge enhancement 
of a computer vision. In this model, u(x,t ) represents an improved version of the 
initial gray level u o(x) of a noisy picture. The anisotropic diffusion div( 1+ |^p ) is 
forward parabolic in the subcritical region where \Du\ < 1 and backward parabolic 
in the supercritical region where \Du\ > 1. 

The expectation of the Perona-Malik model is that disturbances with small gradient 
in the subcritical region will be smoothed out by the forward parabolic diffusion, 
while sharp edges corresponding to large gradient in the supercritical region will 
be enhanced by the backward parabolic equation. Such expected phenomenology has 
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been implemented and observed in some numerical experiments, showing the stability 
and effectiveness of the model. On the other hand, many analytical works have 
shown that the model is highly ill-posed when the initial datum uq is transcritical 
in 0; namely, there are subregions in 0 where \Duq\ < 1 and where |Z>Uo| > 1, 
respectively. For transcritical initial data, due to the backward parabolicity, even a 
proper notion and the existence of well-posed solutions to (11. ip have remained largely 
unsettled. Most analytical works have focused on the study of singular perturbations, 
Young measure solutions, numerical scheme analyses, and examples and properties 
of certain classical solutions; see, e.g., I2i ej m Ha m, nn m\. 

The present paper addresses the analytical issue concerning the existence of certain 
exact weak solutions to problem (11,11) . Let Q? — x (0, T). We say that a Lipschitz 
function u G W 1 '°°(Q,t) is a weak solution to (11.11) provided for all £ G C°°(Ct T ) and 
s G [0,T], 

(—u( t + cr(Du)-D()dxdt — / uo(x)£(x, 0)dx, 

Jn 

where a(p) = (p G R n ) is the Perona-Malik function. The first existence re¬ 

sult on such weak solutions was established by K. Zhang [32j for the one-dimensional 
problem, whose pivotal idea is to reformulate the one-dimensional Perona-Malik equa¬ 
tion as a differential inclusion with linear constraint and then prove the existence using 
a modified method of convex integration following the ideas of [20] [23]. Based on a 
similar approach of differential inclusion, we have recently proved in [19] that for all 
dimensions n if the domain hi is a ball and the nonconstant initial function uq is 
smooth and radially symmetric then (11.11) admits infinitely many radially symmetric 
Lipschitz weak solutions. 

The main purpose of this paper is to extend the results of [19] [32] to problem (11.11) 
on all n-dimensional smooth convex domains for all nonconstant smooth initial data. 

Our main result of the paper is the following theorem. 


(1.3) / u(x, s)£(x, s)dx + 



Theorem 1.1. Let f 2 C R n be a bounded convex domain with dLl of C 2+a and 
let un G C 2+a (PT\ satisfy (|1.2|) for some constant 0 < a < 1. Then in possesses 
infinitely many weak solutions. Moreover, if ||-D«o||l°°(Q) > 1 and A > 0, then these 
weak solutions u will satisfy the almost gradient maximum principle: 

< \\Du 0 \\ L ~ m + A. 

This theorem asserts that the Perona-Malik problem (11.11) admits infinitely many 
Lipschitz weak solutions no matter whether the initial datum is subcritical, super¬ 
critical, or transcritical. 

Existence of classical solutions to Problem (II.ip depends heavily on the initial data 
u 0 . Kawohl & Kutev tra showed that a classical solution exists in any dimension 
if uq is subcritical in O (see also [TS]). Later, Gobbino [15] showed that the problem 
cannot admit a global classical solution when n — 1 if Uq is transcritical. Recently, 
Ghisi & Gobbino [T3J 33] have studied the existence and properties of certain clas¬ 
sical solutions of the Perona-Malik equation in the one-dimensional or n-dimensional 
radially symmetric cases with suitably chosen initial data; their initial values can be 
arbitrarily given in the subcritical region, but the values in the supercritical region 
must be predetermined by the subcritical initial values. 
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We remark that the convexity of the domain is needed to guarantee a gradient max¬ 
imum principle for the classical solution to initial-Neumann boundary value problem 
of a class of quasilinear uniformly parabolic equations (see Theorem 12. II below). This 
gradient maximum principle turns out to be crucial for the proof of main theorem, 
and an example in [1] Theorem 4.1] showed that such a gradient maximum princi¬ 
ple may fail even for heat equation without the convexity of the domain. However, 
domain convexity seemed to be overlooked in m Theorem 6.1]. 

For the proof of Theorem 11.11 in what follows, we assume the initial function uq 
satisfies 


(1.4) 


uo(x)dx = 0, 


since otherwise one can solve solution u of (II. ip with new initial datum uq = uq — 
■pi f n u 0 dx ; then u = u + j^j f Q u 0 dx will solve (II. ip . 

Our proof is based on a crucial generalization of the ideas of [T9] |32 : , [33]. Let us 
discuss this generalization in some details because it exhibits several different features 
from the one-dimensional setup. 

Assume u G 1F 1,00 (H T ) is a weak solution to (II. ip and suppose there exists a vector 
function v G W l,00 {flx] R") such that div v = u and v t = a {Du) a.e. in fl T . Let 
w = (u, v ): Qt R 1+n , with space-time Jacobian matrix denoted by 

'Du uf 
Dv v t 


Vw = 


as an element in the matrix space M^ 1+ ^x( n+1 ). Given s G R, define the set K(s) in 

]yj(l+n)x(n+l) ky 


(1.5) 


K(s) = 


p c 

B a{p) 


p G R", c G R, B G M nxn , tr B = s 


Then w = (u, v ) solves the nonhomogeneous partial differential inclusion: 

X7w(x,t) G K(u(x,t )), a.e. (x,t) G fix- 

Conversely, suppose we have found a function $ = (u*,v*), where u* G kF 1,00 (fLr) 
and v* G W 1,00 (flx] R n ), such that 

u*{x, 0) = Uo(x) (x E ft), 
divv* = u* a.e. in fix, 
v*{-,t) ■ n| 9c = 0 Vt G [0, T], 

Assume w = (u,v) G JF 1 ’ 00 (H'r; R 1+n ) solves the Dirichlet problem of nonhomoge¬ 
neous differential inclusion: 

Vw{x,t ) G K(u(x,t)), a.e. (x,t) G fix, 
w(x,t ) = <&(x, t), (x,t) G dflx■ 


( 1 . 6 ) 


(1.7) 


Then it can be verihed that u is a weak solution to (II. ip (see Lemma 13.21) . 

The Dirichlet problem (II.7K falls into the framework of general nonhomogeneous 
partial differential inclusions studied by Dacorogna & Marcellini cm using 
Ba.ire’s category method and by Muller & Sychev [2S] using the convex inte¬ 
gration method; see also [20] . Study of such differential inclusions has stemmed 
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from the successful understanding of homogeneous differential inclusions of the form 
Du(x) G K first encountered in the study of crystal microstructure by Ball & 
James [2], Chipot & Kinderlehrer [7] and Muller & Sverak plj. Recently, 
the method of differential inclusions has been successfully applied to other important 
problems ; see, e.g., P UTl 12311251 1301 I3l] . 

We point out that the existence result of [26] is not applicable to problem (jl . 7[) 
even in dimension n — 1, as has already been noticed in [32] 03]. A key condition 
in the main existence theorem of [25] . when applied to (ll.Tjl . would require that the 
boundary function <f> satisfy 

V<h(a;, t) G U(u*{x , t)) U K(u*(x , t)), a.e. (x, t) G h2 T , 

where U(s) C M' 1+n ) x b +1 ) (s G R) are bounded sets that are reducible to K(s) in 
the sense that, for every s 0 G R, £ 0 G U(s 0 ), e > 0, and bounded Lipschitz domain 
G C R n+1 , there exist a piecewise affine function w G Wq ,00 (G\ R 1+n ) and a 8 > 0 
satisfying, for a.e. z = (x, t ) G G, 

£o + Vu;(;s)G Q U(s), f dist(£ 0 + Vw(z), K(s 0 )) dz < e\G\. 

|s—s 0 |<<5 


The second condition would imply tr£> 0 = s o for each £o = G U(s 0 ) and 

so G R; but then n| s _ So | <( st/(s) = 0, which makes the first condition impossible. 

However, certain geometric structures of the set K( 0) turn out still useful, especially 
when it comes to the relaxation of homogeneous differential inclusion Vca(z) G K( 0) 
with uj = We explicitly compute the first-order lamination set L(K( 0)) of 

K( 0) consisting of all £ G M (1+n ) x ( n+1 ) \ K( 0) such that £ = A£i + (1 — A)£ 2 for some 
A G (0,1) and £i,£ 2 G K( 0) with rank(£i — £ 2 ) = 1. We obtain the explicit formula 
(see Theorem ED 


L(K( 0)) = 


P 

B 


trR = 0, |£| 2 + (p-(3) 2 -p-p < 0 


which enables us to extract enough information on the diagonal components of dif¬ 
ferential inclusion Vco(z) G AT(0) and establish a relaxation result on (. Dip,i/; t ) (see 
Theorem 14.6p . Although for such relaxation we must have div0 = 0, the resulting p t 
can be arbitrarily small; this is important for the subsequent handling of the linear 
constraint div v — u in problem (II .7|) . 

Another difficulty concerning problem (II .7]) is that when n — 1, one can control 
||i>x||L°°(n T ) in terms of u = v x (see [32]): however, for n > 2, it is impossible to 
control \\Dv\\l°°(q t ) in terms of u — div v. So, if n > 2, the space W 1)00 (yiT] R n ) is 
not suitable for the function v. It turns out that a suitable space for v is the space 
W 1,2 ((0, T); L 2 (h2; R n )) of abstract functions (see Lemma 13.R ; in this setting, the 
linear constraint div v — u must be understood in the sense of distributions. 

We design a new approach to overcome the lack of control on Dv: instead of 
defining an admissible class for w = (u,v), we define a suitable admissible class for 
only the functions u G iy 1 ’ 00 (h2'r); treating v as auxiliary functions. Of course, during 
all the relevant constructions, the linear constraint div v = u must be satisfied. In 
this regard, we need a linear operator 1Z that serves as a (distributional) right inverse 
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of the divergence operator: divAl = Id. By the results of ||4], such an operator may 
not exist as a bounded operator on certain spaces, but for our purpose, it suffices to 
construct such an operator TZ that is bounded from L°°(Q x I) to L°°(Q x /; R n ) for 
the box domains Q x / in R n+1 ; this is achieved by following some construction in 

m- 

Finally we remark that although the result of this paper heavily relies on the explicit 
formula of L(K( 0)), the method can handle some general forward-backward parabolic 
equations; however, we do not intend to discuss further results of this direction in the 
present paper. 

The rest of the paper is organized as follows. In Section 2, we collect several nec¬ 
essary preliminary results, some of which cannot be found in the standard references. 
In Section 3, we set up a new general procedure for proving Theorem 11.11 under a 
pivotal density hypothesis of an admissible class li\ this setup is suitable for a Baire’s 
category method and simplifies some of the arguments even for the one-dimensional 
problem. In Section 4, as the heart of the matter for fulfilling the density hypothe¬ 
sis and thus proving Theorem 11.11 we present the essential geometric considerations, 
including an explicit computation of the set L(K( 0)) above and establishing a crit¬ 
ical relaxation property (Theorem 14.6[) by convex integration with linear constraint. 
In Section 5, we construct the suitable admissible class U after defining a specific 
boundary function $ = (u*,v*). In Section 6, we fulfill the key density hypothesis 
for admissible class U (Theorem 16.111 and finally complete the proof of Theorem 11.11 
according to the setup of Section 3. 

2. Some preliminary results 

2.1. Uniformly parabolic quasilinear equations. We refer to the standard ref¬ 
erences (e.g., |2H|22]) for general theory of parabolic equations, including some no¬ 
tation concerning functions and domains of class C k+a for integer k > 0 and number 
0 < a < 1 . 

Assume / e C 3 ([0, oo)) is a function satisfying 

(2.1) 9 < f(s) + 2sf(s) < 0 Vs > 0, 

where 0 > 9 > 0 are constants. This condition is equivalent to 9 < (s/(s 2 ))' < 0 for 
all s G R; hence, 9 < f(s) <0 for all s > 0. Let 

Mp) = f(\p\ 2 )p ip e R n )- 

Then we have 

A l Pj {p) = fi\p\ 2 )S ij + 2f(\p\ 2 )p i p j (i,j = 1,2, - ■ • ,n; p e R n ) 
and hence the uniform ellipticity condition : 

n 

(2.2) 9\q\ 2 < ^ A ' Pj (P)mj < ©I# Vp, q G R n . 

i,j= 1 

Theorem 2.1. Let fl C R n be a bounded convex domain with dfl of C 2+a and 
uq G C ,2+Q (f2) satisfy Duq • n = 0 on dfl. Then the initial-Neumann boundary value 
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problem 


(2.3) 


{ Ut = di y(A(Du)) 
du/d n = 0 
u(x, 0) = u 0 (x) 


in Dt, 

on dfl x (0, T), 
for x G hi 


has a unique solution u G C 2+Q ’ 2 2 ° (D, T ) ■ Moreover, the gradient maximum principle 
holds: 


(2.4) \\Du\\L°°(n T ) — ||-DmoIU°°(Q). 

Proof. 1. As problem (12. 3p is uniformly parabolic by (12.2[) , the existence of unique 
classical solution u in follows from the standard theory; see [22* The¬ 

orem 13.24], To prove the gradient maximum principle (I2.4H . note that, since A G 
(7 3 (R n ), a standard bootstrap argument based on the regularity theory of linear para¬ 
bolic equations [2T1I22] shows that the solution u has all continuous partial derivatives 
u Xi xjx k and u Xi t within f It for 1 < i, j, k < n. 

2. Let v = |/>u| 2 . Then, within fl T , we compute 

Av = 2 Du ■ D(Au) + 2\D 2 u\ 2 , 
u t = di v(A(Du)) = di v(f(v)Du) = f'{y)Dv ■ Du + f(v)Au, 

Du t =f"(y)(Dv ■ Du)Dv + f\v){D 2 u)Dv 

+ f'(v)(D 2 v)Du + f'(v){Au)Dv + f{v)D(Au). 

Plugging these equations into v t = 2 Du ■ Du t , we obtain 

(2.5) v t - jC(v) - B ■ Dv =-2f(\Du\ 2 )\D 2 u\ 2 <0 in fi T , 
where operator C(v) and coefficient B are defined by 

C(v) = f(\Du\ 2 )Av + 2f'(\Du\ 2 )Du ■ ( D 2 v)Du , 

B = 2f"(v)(Dv ■ Du)Du + 2 f\v)(D 2 u)Du + 2f{v){Au)Du. 

We write C(v) = Y^j =i a ij v x iXj , with coefficients cqj = aij(x,t ) given by 

a ij = Al Pj ( Du ) = f(\ Du \ 2 )5ij + ‘2f\\Du\ 2 )u Xi u Xj ( i,j = 1, • • • ,n). 

Note that on Qt all eigenvalues of the matrix (a tJ ) lie in [0,0], 

3. We show 

max v(x,t) = ma_xt;(x,0), 

(x,t)(EfiT 

which proves (12.40 . We prove this by contradiction. Suppose 

(2.6) M max v(x, t) > maxw(i, 0). 

(x,t)£flT x££l 

Let (x 0 ,t 0 ) G Cl T be such that v(x 0 ,t 0 ) = M; then t 0 > 0. If x 0 G hi, then the strong 
maximum principle applied to (12.5[1 would imply that v is constant on Q <0 , which 
yields v(x, 0) = M on T2, a contradiction to (12.6[) . Consequently xq G dD and thus 
v(xo,t 0 ) = M > v(x,t) for all (x,t) G Or- We can then apply Hopf’s Lemma for 
parabolic equations [28] to (12. 5p to deduce dv(xo,t 0 )/dn > 0. However, a result of 
[H Lemma 2.1] (see also [16.i Theorem 2]) asserts that dv/dn < 0 on dfl x [0, T] 
(convexity of H is used and necessary here), which gives a desired contradiction. □ 
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2.2. Modification of the Perona-Malik function. We need to modify the Perona- 
Malik function a(p) = 1+ P \ p \2 to obtain a uniformly parabolic problem of type (12.31) . 
For this purpose, let 

p(s) = Y^ ( s ^°) 

and, for 0 < 5 < 1/2, let m = m±(6) be the solutions of p(m) = 5; that is, 


(2.7) 


m±(5) 


1 ± y/1 - 4 S 2 
25 


The following result can be proved in a similar way as in [6j [32]; we omit the proof 
(see Figure 1). 


Lemma 2.2. Let 0 < 5 < 1/2 and 1 < A < m + (5). Then there exists a function 
p* G C^QO, oo)) satisfying that 

p*(s) = p(s ) V 0 < s < m-(6), 
p*(s) < p(s) V m-(5) < s < A, 

6 < (p*)'(s) <0 V0<s<oo 

for some constants 0 > 6 > 0. Moreover, define /(0) = 1 and f(s ) = p*(\/s)/\fs for 
s > 0; then f G C 3 ([0, oo)) and (12.ip is fulfilled. 



FIGURE 1. The graphs of function p(s) and the modihed function p*(s) 
in Lemma [2.21 

2.3. Right inverse of the divergence operator. To deal with the linear constraint 
div v = u, we follow an argument of [4j Lemma 4] to construct a right inverse 1Z of 
the divergence operator: div 1Z = Id (in the sense of distributions in fir)- For the 
purpose of this paper, the construction of TZ is restricted to the box domains, by which 
we mean domains given by Q = J\ x J 2 x • • • x J n , where Jj = (a;, bf) C R is a finite 
open interval. 

Given such a box Q, we define a linear operator 7 Z n : L°°(Q ) -A L°°(Q ; R”) induc¬ 
tively on dimension n. If n = 1, for u G L°°(Ji), we define v = IZ\U by 

rx i 

v(xi) = / u(s)ds (xi G Ji). 

J ai 
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Assume n = 2. Let u G L°°(J 1 x J 2 ). Set u(x i) = f^u(xi,s)ds for Xi G Ji. Then 
u G L°°(Ji). Let b = IZiU] that is, 

j-x 1 /-Xi r-b 2 

b(xi) = / u(s)ds = / / u(s,T)drds (xi G Ji). 

J a\ J a\ J CL2 

Let p 2 ^ Cf°(a 2 ,b 2 ) be such that 0 < p 2 (s) < ^Q— and p 2 (s)ds = 1. Define 
r? = 7^ 2 u G L°°(Ji x J 2 ;R 2 ) by u = (u/r; 2 ) with v 1 (xi,x 2 ) = p 2 (x 2 )v(xi) and 

rx2 rx2 

v 2 (x\, x 2 ) = / u(xi, s)ds — u(xi) / p 2 (s)ds. 

J CL2 J CL2 

Note that if u G bb 1,0O (Ji x J 2 ) then u G bb 1,0O (Ji); hence u = lZ 2 u G hL 1,00 (Ji x 
J 2 ; R 2 ) and divu = u a.e. in ,/ x J 2 . Moreover, if u G C 1 (Ji x J 2 ) then u is in 

c\l^r 2 - r 2 ). 

Assume that we have defined the operator 1Z n -\. Let u G L°°(Q ) with Q = 
Ji x J 2 x • • • x J n and x = ( xx n ) G Q, where x' G Q' — J\ x ■ ■ • x J n _! and x n G J n . 
Set u(x') = j^ n u(x',s)ds for x' G Q’. Then u G L°°(Q'). By the assumption, v = 
Tl n -iu G R n_1 ) is dehned. Write v(x') = (Z 1 (x l ), • • • , Z n_1 (x')), and let p n G 

Cf°(a n , b n ) be a function satisfying 0 < p n (s ) < fc c ° and J^ n p n (s)ds = 1. Dehne v = 
lZ n u G L°°(Q ; R") as follows. For x = (x',x n ) G Q, u(x) = (u 1 (x), u 2 (x), • • • ,u n (x)) 
is dehned by 

u fc (x', x n ) = p n (x n )Z k (x') (k = 1, 2, ■ ■ • , n - 1), 

PXn PX n 

v n (x', x n )= / u{x', s)ds — u{x') / p n (s)ds. 

J d n J CL n 

Then 7T„: L°°(Q) —>■ L°°(Q ; R n ) is a well-dehned linear operator; moreover, 

(2.8) ||7?. n 'u||£°o(Q) < C n (|</11 + • • • + | Ail)|| w||l°°(q), 

where C n > 0 is a constant depending only on n. 

As in the case n = 2, we see that if u G hF 1,00 ((5) then v = 7 Z n u G TW 1,00 (Q; R n ) 
and divu = u a.e. in Q. Also, if u G C l (Q) then v = TZ„u is in C 2 (Q; R n ). Moreover, 
if u G Wi'°°(Q) satishes f Q u(x)dx = 0, then one can easily show that v = TZ n u G 

WS’°°(Q-, R"). 

Let / be a finite open interval in R. We now extend the operator 7 Z n to an operator 
1Z on L°°(Q x I) by defining, for a.e. (x,t) G Q x /, 

(2.9) (Jlu)(x,t) = (Jl n u(-,t))(x) V uG T°°(Q x /). 

Then 7?.: L°°(<3 x I) —» L°°(Q x /; R") is a bounded linear operator. 

We have the following result. 

Theorem 2.3. Let u G hF 0 1,oo ((5 x /) satisfy j Q u(x,t) dx = 0 /or all t G /. Then 
u = 7Tu G W 0 1 ’ oo (Q x /; R n ), div v = u a.e. in Q x I, and 

(2.10) ||^||l°°(Qx/) < Cn (|^l| + • • • + |T n |)||u t ||£oo(Q x j), 

where Q = J\ x • ■ ■ x J n and C n is the same constant as in (12.8[) . Moreover, if 
u G C 1 (Q x I) then v = TZu G C l (Q x /; R n ). 
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Proof. Given u G Wq’°°(Q x /), let v — 7 Zu. We easily verify that v is Lipschitz 
continuous in t and hence v t exists. It also follows that v t = 7 Z(u t ). Clearly, if 
fgu(x,t)dx = 0 then v(x,t) = 0 whenever t E dl or x G dQ. This proves 

v G W^°°(Q x /; R ra ) and the estimate (12. 10p follows from (12.811 . Finally, from the 
definition of 7 Zu, we see that if u G CT{Q x I) then v = 7 Zu eC 1 (Q x I; R n ). □ 


3. General setup for existence 

In this section we set up the general procedure for proving our main theorem, 
Theorem 11.11 


3.1. Sufficient conditions for weak solutions. Since our setup differs from the 
usual formulation of differential inclusions, we first prove the next two results to 
clarify some relevant issues, which are elementary but not too obvious. 

Lemma 3.1. Supposeu G lF 1 ’ 00 (f2r) is such thatu(x, 0) = uq(x) (x G fl), there exists 
a vector function v G 1T 1,2 ((0, T); L 2 {Q; R n )) with weak time-derivative v t satisfying 
v t = a {Du) a.e. in VLt, o,nd for each ( G C°°(f2 r ) and t G [0, T ], 


(3.1) / v(x, t) ■ D({x, t) dx = — / u(x, t)((x, t) dx. 

Jn Jn 

Then u is a weak solution to (11.11) . 

Proof. To verify (11.31) . given any ( G C°°{Cl T ), let 

g(t) — / u(x,t)C(x,t)dx, h{t) — / u(x,t)( t (x,t)dx (0 <t<T). 


Then for each if G C“(0, T), by (13.11) . 


ift{t)g{t) dt — — 


if{t)h{t ) dt = 


n 

>o Jn 

fj 

Jo J fi 


if t (t)v(x,t) ■ D({x,t ) dxdt , 
if(t)v(x, t ) • D( t {x, t ) dxdt. 


Since v G 1T 1,2 ((0, T); L 2 (f2; R”)) and v t = a{Du), one has 


(■ ip{t)D((x , t)) t ■ v(x, t)dxdt = 


if{t)a(Du(x, t )) • D({x , t) dxdt. 


'o Jn 


'o Jn 


Now as {ifD(f) t = if t D( + ipD(f t , combining previous equations, we have 

[ if t (t)g(t) dt — f if{t)(—h(t)+ f a{Du(x,t)) ■ D((x,t) dx\dt; 
Jo Jo \ Jo. J 

this proves that g is weakly differentiable in (0, T ) with weak derivative 


g\t) = li(t) — / a{Du{x,t)) ■ DC,{x,t) dx a.e. i G (0, T). 
Jn 

From this, upon integrating, (II.3ft follows for all s G [0,T]. 


□ 
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Let $ = («*, v*) G W 1)CO (VL t \ R 1+n ) satisfy (US]), and let W/l°°(n T ) and Wy*°°(0 T : R n ) 
denote the usual Dirichlet classes with boundary traces u*, v*, respectively. 

Let U be some nonempty and bounded subset of W^°°(Rr) such that for each 
u 6 if, there exists a vector function v G R n ) satisfying divu = u a.e. in 

fbr and ||ut||L°°(n T ) <1/2. Any such set U is called an admissible class. 

Given any e > 0, define U e to be the set of u G U such that there exists a vector 
function v G R") satisfying divv = u a.e. in Qt, ||'n’t||z,°°(n r ) < 1/2, and 


I v t (x,t) 


cr(Du(x,t))\ dxdt < e\D t 


J Qj 1 

(Note that U t =U for all e > 1.) 

Lemma 3.2. Let u G U.. Then any vector function v G Wj*°°(Rr; R n ) determined 
above satisfies the integral identity (13.II) for each Q G C°°{VLt) and t G [0 ,T\. 

Proof. Let ( G C°°(D T ) and define 


h(t) — / (v(x, t) ■ Df(x,t) + u(x, t)((x,t)) dx. 


Then h is continuous on [0, T] and for each if G C 1 [0, T], 


h(t)ip(t) dt = 


fs 

n 

'o Jn 

ff 

r-T r 
1 0 JCl 


i/j(t) ( v(x , t ) • D((x, t) + u(x, t)((x, t )) dxdt 
ijj(t) (v(x, t ) • D((x, t ) + divu(a;, t)((x, t )) dxdt 
cliv(</(;r, t)y(t)v(x, t )) dxdt 
div(C(x, t)'if{t)v*(x 1 1 )) dxdt = 0, 


resulting from v\qq t = v*\qci t and v*(-,t) ■ n|gn = 0 for all t G [0, T\. ffence h = 0 on 
[0,T]. This completes the proof. □ 

3.2. General existence theorem by Baire’s category method. We prove a 
general existence theorem under a density hypothesis. 


Theorem 3.3. LetU C W'/i°°(Ot) be an admissible class. Assume, for each e > 0, 
K is dense in U under the L°°-norm. Then, given any ip G hi, for each p > 0, there 
exists a weak solution u G to problem (11.11) satisfying ||u — </ 2 ||l°°(q t ) < g. 

Furthermore, if U contains a function in lb/* 00 (Or) that is not a weak solution to 
(ED, then (11.11) admits infinitely many weak solutions. 


Proof. 1. Let X be the closure of U in the metric space L°°(Rr). Then ( X,L °°) is a 
complete metric space. By assumption, U e is a dense subset of X. Furthermore, since 
U is bounded in bb/l°°(fLr), we have X C W/i°°(fbr)- 

2. Let y = L 1 (h2r;R ri ). For h > 0, define T/: X —)■ y as follows. Given any 
u G A", write u = u* + w with w G bF 0 1 ’ oo (r2 r ) and define 

T h (u) = Du* + D(p h *w), 











CONVEX INTEGRATION AND THE PERONA-MALIK EQUATION 


11 


where Ph(z) = h~ N p(z/h), with 2 = (x,t) and N — n + 1, is the standard modifier on 
R iV , and ph * w is the usual convolution on R A with w extended to be zero outside 
VLt- Then, for each h > 0, the map 7) t : (X, L °°) —» ((V, L 1 ) is continuous, and for each 
u E X, 


lim \\T h (u) - Du\\ L irn T ) = lim \\p h * Dw - Dw ||z,i(n r ) = 0. 

h— >0+ h— >-0+ 

Therefore, the spatial gradient operator D : X —> y is the pointwise limit of a sequence 
of continuous functions T/c X -A hence D : X —> y is a Baire-one function. By 
Baire’s category theorem (e.g., [5] Theorem 10.13]), there exists a residual set Q G X 
such that the operator D is continuous at each point of Q. Since X \ Q is of first 
category , the set Q is dense in X. Therefore, given any tp E X, for each p > 0, there 
exists a function u E Q such that ||w — <p\\l° °(n T ) < P- 

3. We now prove that each function u E Q is a weak solution to (11.11) . Let u E Q 
be given. By the density of U € in ( X,L°° ), for each j E N, there exists a function 
Uj E U\/j such that || Uj—u || L oo(q t ) < l/j. Since the operator D: (X, L°°) -E fiy^L 1 ) is 
continuous at u, we have Duj -A Du in R"). Furthermore, from the definition 

of Ui/j, there exists a vector function Vj E R n ) such that, for each ( E 

and t E [0, T], 


Vj(x,t) ■ D((x,t)dx — — / Uj(xfi)((x,t) dx, 


(3.2) 


(^j)tII t°°(Ut) < 2’ 


'o Jn 


| (vj)t — cr(Duj) | dxdt < -|fhr| 


Since Vj(x, 0) = v*(x, 0) E W l, °°{yi\ R n ), from ||(t'y)«||< 1/2, it follows that 
both sequences {vj} and {( 13 )*} are bounded in L 2 (VL T ] R n ) ~ L 2 ((0, T); L 2 (h2; R n )). 
We may assume Vj — ^ v and (vj) t — " v t weakly in L 2 ((0, T); L 2 (h2; R n )) for some 
v E 1F 1 , 2 ((0, T); L 2 (D] R n )). Upon taking the limit as j —> 00 in (13. 2 p and noticing 
v E C'([0, T]; L 2 (h2; R n )), we obtain that 


/ v(x, t)-D((x, t) dx — — / u(x, t)C(x, t) dx (f 6 [0,T]), 

Jn Jn 

Vt(x,t) = a(Du(x,t )) a.e. (x,t) E Dt- 

Consequently, by Lemma 13.11 u is a weak solution to (II.ip . 

4. Finally, assume U contains a function that is not a weak solution to (II. 1 [) : hence 
Q U. Then Q cannot be a finite set since otherwise the L°°-closure X = Q = U 
would be a finite set, making U = Q\ therefore, in this case, ( 11 . 11 ) admits infinitely 
many weak solutions. This completes the proof. □ 

The rest of the paper is devoted to the construction of a suitable admissible class 
U C PF 1 , 00 (fFr) fulfilling the density property: 


(3.3) 


lA t is dense in U under the L°°-norm for each e > 0. 
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4. Geometric considerations: Relaxation of Vlu(z) g K( 0) 

Let AT(s) be the matrix set defined by (1 1.5 j) above. Since AT(s) is a translation of 
set K( 0), we focus on the set K 0 = A'(0); that is, 


AV, = 


P 


B <j{p) 


p g R", c g R, B G M IlXn , trR = 0 


where a(p) = y+j^p is the Perona-Malik function. 


4.1. Rank-one lamination of A' (l . We first compute certain rank-one structures of 
the set A'q. 

Let A(A' 0 ) be the set of all matrices £ G Md +n ) x ( n+1 ) that are not in K 0 but are 
representable by £ = A£i + (1 — A)£ 2 for some A G (0,1) and £i, £ 2 G K 0 with 
rank(£! — £ 2 ) = 1 , or equivalently, 

A(A'o) = {£ ^ A'q | £ + t±rj G A ' 0 for some t_ < 0 < £+ and rank?/ = 1}. 

Suppose £ = ^ ^ G A(A'o), with £ + t±rj G A' 0 , where t_ < 0 < £+ and p is a 
rank-one matrix given by 

= (“) ® M) = ( a ®, it) ' ° 2 + |o l 2 ^ °- 62 + l"l 2 ^ °- 

for some a, b G R and a, q E R n ; here a ® q denotes the rank-one or zero matrix 
(ot l qj) in M nxri . 

Condition £ + t±p G K 0 with f_ < 0 < t + is equivalent to the following: 

(4.1) trR = 0, a ■ q — 0, cr(p + t±aq) = (3 + t±ba. 

If aq = 0, then a(p) = (3 + tba has two different solutions of t only when ba = 0, but 
then we would have a(p) = (3 and thus £ G A' 0 , a contradiction. Therefore, aq ^ 0. 
By rescaling r/ and t±, we assume a = 1 and \q\ = 1; namely, 

t)- = = 


Case 1. Assume ba = 0. In this case, by (14. 1[) , the equation a(p + tq) = {3 has two 
solutions of t of opposite signs and thus we must have p — xq and /3 = uq, and 
cr(xq + tq) = uq becomes a quadratic equation x + t = u + u(x 2 + 2 xt + t 2 ), which 
has two solutions t — t± of opposite signs if and only if u 7 ^ 0 and x 2 — ^ + 1 < 0 ; 
this condition can be written as 


\/3\ 2 + (p ■ (3) 2 — p ■ f3 = u 2 + (xu) 2 — xu < 0 . 


Remark 4.1. In this case one can always select rj 


q b 
0 0 


1, b G R. This is 


the case for the one-dimensional problems studied in [19] [32], :33], where the existence 
results are primarily proved based on the structure of such 77 ’s. However, if n > 2, 
such 7 /’s are not sufficient to characterize all the rank-one structures. Case 2 below 
thus becomes pivotal. 
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Case 2. Assume ba ^ 0; so b ^ 0 and a/0. In this case, we write 


’ ,= (b®«%)’ = '-9 = °, 7#°, M°- 

Since the equation a(p + tq) = P + i 7 has two solutions t — t± of opposite signs, we 
must have p = xq + y'y and (3 = uq + V'y, and the equation cr(p + tq) = P + t'y becomes 
a system of two equations: 


x + t = u{l + (x + t) 2 + | 7 | 2 y 2 ), 

y = (u+f)(i + (*+f ) 2 + | 7 iy). 


This system has two solutions t — t± of opposite signs, and thus and y ^ 0. So 

(14.2p is equivalent to a system of two quadratic equations: 


(4.3) 


t 2 + ( 2 x - i)t + x 2 + | 7 |V + 1-2 = 0 , 
f 2 + (u + x)i + xu — yu = 0 . 


The necessary and sufficient condition for (14.3p to have two solutions t = t± of 
opposite signs is that the two quadratic equations of t have the same coefficients and 
the constant terms are negative, which yields that 

1 o i 12 2 -i X r\ 

x = —hr, x + 7 j/" + l- = xv — yu < 0 . 

u u 

Here, if v — 0, then x — and taking this into the inequality, we have 1 + | 7 | 2 y 2 < 0, 
a contradiction. So v 0. Therefore 

/ . . \ XU 1 

(4.4) uv = xu — 1, | 7 | = 


yv 


r 


and 

(4.5) 


XV 


yu =(£-£) {XU - 1 ) < 0 . 


We now solve for x, y, u, v from (14.4H in terms of p and p. From p = xq + yj, (3 = 
uq + vj, it follows that 


(4.6) q = 

xv — yu 

By (14.411 and (14.611 . we have 


-( vp-y/3 ), 7 = 


xv — yu 


-up + x(3). 


where • (3 
(4.7) 

Moreover, 


p ■ /3 = xu + yv\^f\ 2 


1 


V X 

y u 


l\p\ 2 -p-P 

p-p-\pr 


\P\ 2 7 ^ 0 by (14.511 . Let k = x/u, l = y/v. Then 


1 = (1 ~ P ■ /3)\p\ 2 - p ■ (3 

1-P-/+ (l-p-p)p-p-\P\ 2 ' 


x _ y _ \p\ 2 -ip-P _, = \p-W\ 2 

U V p ■ P — l\P\ 2 p ■ P — l\P\ 2 ' 
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From |g| = 1, we have xv — yu = {k — l)uv = —|vp — lv/3\ and hence 

U oA f \\P~W\ , ^P' P~PP I 2 

(4.8) -sgnty)——— = -sgn(t>) ]p _ m , 

where sgn(n) = v/\v\ is the sign of v ^ 0. We then obtain x, v, y by 


(4.9) 


x = kn, v = x-- = kn-~, 
u u 


y — Iv — Iku — 


l 


u 


In this way, we have solved x, y, u, v in terms of p, f3, uniquely up to the sign change. 
We can check that both conditions in (14.41) are satisfied. 

Let us consider inequality (j4.5j) for these solutions. Equation on |q | 2 in (14. 4 p implies 
-(1 — xu ) > 1 and hence the inequality (14. 5 p yields -(1 — xu ) > -(1 — xu) > 1. So 
0 < xu < 1 and thus x/u > y/v > 1, i.e., k > l > 1. Then we deduce the inequality 

(4.10) |/3| 2 + (p - /3) 2 -p - ^ < 0. 

4.2. Exact formula of L(Kq). In fact, inequality (14.101) exactly characterizes the 
set L(Ko). We have the following result. 

Theorem 4.1. 


(4.11) 


L(K„) = 


p c 

B 0 


tr B = 0, \/3\ 2 + (p-/3) 2 -p-p < 0 


Moreover, given any £ G L(Kq), there exist a rank-one matrix 

q h 


V = 


b 


q 7 


with \q\ = 1 , 7 • q — 0 , b ^ 0 and two numbers t- < 0 < t + such that 

£ + t±rj G Ao, 

where | 6 | > 0 can be arbitrarily small. 

Proof. Let S be the set defined on the right-hand side of (14.lip . The previous cal¬ 
culations show that L(K 0 ) C S. To verify the reverse inclusion S C L(K 0 ), let 


£ = 


G S. Then 


+ ('P ' /3 ) 2 ~ P ■ P < 0. So 0 < p ■ P < 1 and 


p c 

B b 

(1 — p ■ P)p ■ P — \P \ 2 > 0, and hence we can define /, k by (14.7j) . so that l > 0, 
k > 0. To fix the sign, we define u by (14.8p with + sign: 


u = 


P-P- 


(1 -P- P)P- P~ I P? 


> 0 . 


\p-ip\ \{l-p-P)p-P\ 

We now define x,v,y by (14.Op . Then x > 0, v < 0, and y < 0. After deducing that 
xv — yu < 0, we finally define q, 7 by (14.61) . It is then straightforward to check the 
following: 

p = xq + yr/, P = uq + v'f, \q\ = 1 , 7 ■ q = 0, 


x — — \-v, 
u 


x 2 + 


l2 y 2 + 1 - — xv 

u 

In particular, equation cr(p + tq ) = P + tj has two solutions t = t± with f_ < 0 < t+. 
q b N 

i7 ®q 7, 

£ G L(K 0 ). The proof is now complete. 


x 


yu < 0 . 


Now let r) = 


, where b 7 ^ 0 is arbitrary. Then £ + t±rj G K 0 , and so 

□ 
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Remark 4.2. The quantities q G S”” 1 , 7 G R ra and t± defined in the proof depend 
continuously on (p, (3) G R n+ " with |/3| 2 + (p ■ /3 ) 2 -p-/3< 0. We may also take b > 0 
(or b < 0 ) to be a continuous function of all such (p,/?). 


4.3. The approximating sets Ss and S 5 . Given any 0 < <5 < 1/2, let 
(4.12) = {(p, /3) G R n+n | 51 (1 - p ■ P)p - p\ + \P\ 2 + (p • P) 2 - p • p < 0}, 


5 ,= 


p c 
B P 


trR = 0, (p,/3) G S 5 


Then S'o = L(K 0 ). Immediate properties of the open sets Ss are that 
Ss 2 C Ss 1 for 0 < hi < S 2 < 1/2, S 1/2 = 0, 

S s ^<b for 0 < S < 1 / 2 , S 0 = (J 5 5 . 

0<5<l/2 

In what follows we always assume 0 < S < 1/2 unless otherwise stated. 


Proposition 4.2. Let £ G Ss and £± G K 0 with rank(£ + — £_) = 1 satisfy that £ lies 
in the open line segment (£_,£+). Then (£_,£+) C Ss- 


Proof. Consider functions 

f(s) —1(1—p■ @)p — p\, G(e) = i/ 3| 2 + fp.« 2 - P -ft C X 

Then both F and G vanish on set Kq. For the given £ and £±, let f(r) = P(£ r ) and 
g{r) = G(£ T ), where £ T = r£ + + (1 — r)£_. The rank-one condition implies that the 
corresponding term p T ■ is linear in r; hence g{r) is a quadratic polynomial of r 
and /(r) is the length of a vector quadratic in r. Since both /(r) and g{r) vanish 
when r = 0 and 1, we must have g{r) = Cir(l — r) and /(r) = G 2 t(1 — r) for some 
constants Ci, C 2 . Since £ = £ A G S'a for some 0 < A < 1, we have 

5/(A) + n(A) = A(1 - A ){SC 2 + Ci) < 0. 

This implies the constant hC 2 + Ci is negative. Hence Sf(r) + g(r) = r(l — r)(hC 2 + 
Ci) < 0 for all r G (0,1), which proves exactly £ T G Ss for all r G (0,1). □ 

Remark 4.3. The lamination convex hull K lc of a set K C M mxn is defined to be 
the smallest set S C M mxn containing K with the property that if £± G S with 
rank(£ + — £_) = 1 then (£_, £ + ) C S (see [9]). As in the proof of Proposition 14.21 one 
can see that Kjf = A'o U L(K 0 ). 

Lemma 4.3. Let m±(5) be defined by (j2.7[) above. It follows that 

sup \p\ = m + (5), inf |p| = 

( p,/3)es s (p,P)£Ss 

Proof. Let (p, P) G Ss- Then 0 < p ■ P < 1 and P ■ ((1 — p ■ P)p — P) > 0. We write 

p- P= \p\\P\ cos#, p- ((l-p-p)p-P) = \P\\(1— p ■ P)p — P\ cosd' 

for some 0 <9, 9' < 7 t/2 . A simple geometry shows that 9 < 9'. As |/3 | 2 + (p • /3 ) 2 — 
p ■ P < 0 , we have 

\/3\ 2 + |p| 2 |/3 | 2 cos 2 6 > - \p\\P\ cos 9 < 0 , 
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and so, 


\p\ < - - — 4 ^ 2 < \p\ < 1 + V 1 4 !gl! 

|P| 2’ 2|/3| cos 9 2|/3| cos 6 ' 


Condition (p,/?) e «S^ becomes 5 < |/3| cos9'] so |/3| > |/3| cos# > 5 and 

l-yT^IF ^ 1- Vl-4|/3| 2 < 1-Vl-41/31 2 


25 


2|/5| 


2|/3| cos6* 


< | P | < EthLASI < 1 + v/rT7iF 


2\/3\ cos 6 


This proves 
(4.13) 


1 - VI - 45 2 


, < inf \p\ < snp |p| < 

45 (p,P)es s (p,/3)es s 45 

Next, fix any (3 e R n with 5 < \f}\ < 1/2. Let 1,1' be any numbers satisfying 

1 


25 

1 + Vl - 45 2 


yi - wv <l<e< i+ yi-tw 


(414) 2|/3| . 2|/J| 

and p = |4/3, V = jV/L Then (p,(3), (p',/3) are both in Sj with \p\ = /, |p'| = l so 

inf Ipl < l < l 1 < sup Ipl. 

(p^V (P,/3)G5 5 

As /, l' are arbitrary and satisfy (I4.14j) . we have 

• t , l-yi-41/jp _ l + v /l-4|^|2 ^ 

, sE IpI £ —Era- < —biTS-£ su p N- 

(p,P)eS s 2|p| 2|p| (p,/3)e5« 

Finally, taking |/3| —)■ 5 + and combining with (j4.13(1 complete the proof. □ 

As an immediate consequence of the previous lemma, we have 


Corollary 4.4. Ss C {(p, ft) | m_(5) < \p\ < m+(5), 5 < |/3| < 1/2}. 

4.4. A useful convex integration lemma. The following result is important for 
convex integration with linear constraint. For a more general case, see [29[ Lemma 
2.1], 

Lemma 4.5. Let Ai, A 2 > 0 and ip = — Aip, r/ 2 = A 2 p with 

’ ,= (h®9 7 )’ 191 = 1 7'« = °. 

Let G C R n+1 5e a bounded domain. Then for each e > 0, there exists a function uj = 
{ip, if) e C c °°(R n+1 ;R 1+ri ) with supp(cj) CC G that satisfies the following properties: 

(a) divV = 0 in G, 

(b) \{z G G | Vuj(z) £ {rji, p 2 »| < e, 

(c) dist(VcW, [771 , 772 ]) < e for all z E G, 
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Proof. 1. The proof follows a simplified version of [29, Lemma 2.1]. Define a map 
V : C' 1 (R n+1 ) —» C°(R n+1 ; R 1+n ) by setting V{h) = (it, i>), where, for h(x,t) G 
C 1 (R n+1 ), 

it (a;, t) = q ■ Dh(x, t), v(x, t ) = — 7 ) Dh(x , f). 

We easily see that R(/i) = (it, v) G C^°(R n+1 ; R 1+n ), suppfRfh)) C supp(h), divu = 
0, and f Rn u(x,t) dx = 0 for all t G R, for all h G C^°(R" +1 ). For h(x,t ) = f{q ■ 
x + bt ) with / G C°°(R), in = (it, i;) = V(h) is given by u(x,t) = f'(q - x + bt) and 
v(x,t) = f'(q ■ x + bt)f, and hence X7w(x,t) = f"(q ■ x + We also note that 
V(gh) = gV(h ) + W?(s) and hence 

(4.15) V?( 5 /i) = gVV{h) + hVV(g) + B{Vg, Wh) V g, h G C°°(R n+1 ), 

where B(yg,Vh) is a bilinear map of Vg and V/t; so |H(V/i, Vg)| < C\X7h\\X7g\ for 
some constant C > 0 . 

2 . Let G e CC G be a smooth sub-domain such that |G \ G e \ < e/ 2 , and let 
p e G Cf°(G) be a cut-off function satisfying 0 < p e < 1 in G, p e — 1 on G e . As G is 
bounded, G C {(x,t) \ k < q- x + bt <1} for some numbers k < l. For each r > 0, we 
can find a function f T G C™(k,l) satisfying 

— Al < Sr < A 2 , \{s G (k,l) | f"(s) { —Ai, A 2 }}| < T, \\fr || L°° + ||/}||l°° < T. 

3. Define u> = = V(p e (x,t)h T (x,t)), where h T (x,t ) = f T (q ■ x + bt). Then 

|| /irHe 1 < G || Sr || c 1 < Gr, uj G G/°(R" +1 ; R 1+n ), supp(cu) C supp(p e ) CC G, and (a) 
and (e) are satisfied. Note that 

M < \p e \\P{h T )\ + \h T \\V(p e )\ < C e T , 

where C e > 0 is a constant depending on 11 p e 11 c ? 1 (<?) - So we can choose a 7 > 0 so 
small that (d) is satisfied for all 0 < r < 77 . Note also that 

{z G G | Vcu(^) ^ { 71 , 72 }} A (G \ G e ) U {z E G e \ f"(q - x + bt) {—Ai, A 2 }}. 

Since \{z G G e \ f” (q ■ x + bt) {—Ai, A 2 }| < N\{s G (k,l) \ f"(s ) </ {—Ai, A 2 }}| for 
some constant N > 0 depending only on set G, there exists a r 2 > 0 such that 

\{z G G | Wuj(z) (£ { 71 , 1 / 2 }}) < | + Nr < e 

for all 0 < r < r 2 . Therefore, (b) is satisfied. Finally, note that 

p e W(hr(x,t)) = pef"(q ■ x + bt )7 G [ 71 , 72 ] in G 

and, by (14.151) . for all z = (. x,t ) G G, 

\Vu(z) - p e VV(h T (x,t))\ < \h T \\VV(p e )\ + \B{Vh T ,Vp e )\ <Cj<e 


for all 0 < r < r 3 , where G' > 0 is a constant depending on ||p e ||c 2 (G), an d r 3 > 0 is 
another constant. Hence (c) is satisfied. Taking 0 < r < min{7T, r 2 , r 3 }, the proof is 
complete. □ 
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4.5. Relaxation of Vc o(z) G K 0 . We now study the relaxation of homogeneous 
inclusion Vui(z) G K 0 . We prove the following result in a form suitable for later use. 

Theorem 4.6. Let 1C be a compact subset of Ss- Let Q x I be a box in R n+1 . Then, 
given any e > 0 , there exists a po > 0 such that for each box Q x / C Q x /, point 
(' p,P ) G 1C, and number p > 0 sufficiently small, there exists a function uj = ((p,i/)) G 
Cf°(Q x /; R 1+n ) satisfying the following properties: 

(a) div -0 = 0 in Q x I, 

(b) (j/ + D(p(z),P' + i>t{z)) G S s for all z G Q x I and \(p',P') - (p, P)\ < po, 

(c) \\u>\\l°°(QxI) < P, 

( d ) Iqxi \P + M z ) ~ a< CP + D<p(z))\dz < e\Q x I\/\Q x J|, 

(e) fg <p(x, t)dx = 0 for all t G /, 

(f) ll^tl|l/°°(Qx/) < P- 


Proof. Let £ = f(p,P) = 


P 

O 


for (p,p) G 1C C Ss, where O is the nxn zero 


matrix. We omit the dependence on (p, P) in the following whenever it is clear from 
the context. Since £ G Ss C L(Kq) on 1C, it follows from Theorem 14.11 and its remark 
that given any p > 0, there exist continuous functions q : 1C —> S n_1 , 7 : 1C —> R n , 
: 1C —> R, and b : 1C —* (0, 00 ) with 7 • q = 0, f_ < 0 < i + on 1C such that letting 
q b N 

~ 9 7, 


V = 


h 


on /C, we have 


£ + t±r) G A' 0 , 0 < b < 


P 


t - 4 




on 1C. 


Writing £± = ^ = £ + t±p on /C, we have £ = A£+ + (1 - A)£_, A = e 

(0,1) on 1C. 

Proposition 0T2] implies that on 1C, both £L = r£++(l—r)£_ and £+ = (1—r)£++r£_ 
belong to Ss for each r G (0,1). Let 0 < r < min^ min {A, 1 — A} < | be a small 
number to be selected later. Let A' = on 1C. Then A' G (0,1), £ = A'£+ + (l —A')£L 
on/C. Moreover, on/C, £^-££ = (l- 2 r)(£ + -£_) is rank-one, [£!,££] C (£_,£+) C S s , 
and 

cr < |£+ - £+| = |£L - £_| = r|£+ - £_| = r(t+ - t_)M < Ct, 

where C = max^(t + — t_)\rj\ > niin^(t + — t_)|? 7 | = c > 0. By continuity, H T = 
U(p,/ 3 )ga:[£L(p, P), C+ip, P)\ is a compact subset of Ss, where Ss is open in 



So d T = c\ist(H T ,d\z 0 Ss) > 0 , where <9|v 0 is the relative boundary in the space £ 0 . 

Let 771 = — Ai 7 = — A'(l - 2 r)(t+ - tJ)p, r) 2 = A 2 t? = (1 - A')(l - 2 r){t+ - tJ)p 
on 1C, where Ai = r(—t+) + (1 — r)(— 1 _) > 0, A 2 = (1 — r)t+ + rt_ > 0 on /C with 
r > 0 sufficiently small. Applying Lemma [4.51 to matrices 771 , r / 2 and set G = Q x I, 
we obtain that for each p > 0, there exists a function uj = (p, if) G Cf°(Q x /; R 1+n ) 
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and an open set G p CC Q x / satisfying the following conditions: 


(4.16) 


ru) 

div-0 = 0 in Q x /, 


(2) 

|(Q x I)\G P \ <p; e + V(n(z; 

G {£1, ££} for all zGG 

(3) 

£ + Vc j(z) G [£I,£+]p for all z 

GQxJ, 

(4) 

1 M L°°(QxJ) < Pi 


(5) 

Jq <p(x, t) dx — 0 for all t G I, 


k (6) 

11 Pt 11 Z.°° (Q x /) < 2P, 



where [£L>£+]p denotes the p-neighborhood of closed line segment [£!,£+]. From 
(14.161 3). (14.161 61 follows as 

\<Pt\ < |C+ - c_| + p = (t+ - t-)\b\ + p < 2p in Q x I. 

By (14.161 3). |/3 + 'ipt(z)\ < C + p for z G Q x J; hence 

/ |/3 + 'ip t - u(p + Dp)\dz 


(4.17) 


J Qxl 


|/3 + — cr(p + D<p)\dz + (C + p + -)p 


< \Q x I\ max{|/3± - a(p±)|} + (C' + p+ -)p 

< C\Q x J|r + |Q x I\ max{|cr(p±) - <x(p±)|} + (C + p + ^)p, 


where f 


± 


Note (a), (c), (e), and (f) follow from (14. 16|) . where 2 p in (14.161 61 can be adjusted to 
p as in (f). By the uniform continuity of cr on the set J = {p' G R n | \p'\ < m + (8) + 1}, 

whenever p', p" G J and 


3|Qx/| 


we can find a p' > 0 such that | cr(p') — cr(p")| < 

|p' — p"| < p', where m + (8) > 0 is the number defined in Lemma 14.31 We then 
choose a r > 0 so small that Cr < p' and C\Q x J|r < |. Since p±, p T ± G J and 
| p± —p T ± | < Cr < p', it follows from ()4.17(1 that (d) holds for any choice of p > 0 with 
(C + p + \)p < ^ Q*j| . Next, we choose a po > 0 such that po < ^. If 0 < p < po, 
then by (14.161 1) and (14.161 3). for all z G Q x I and | (p',/3') — (p, /3)| < p 0 , 

f (p',/3') + X7u{z) G So, dist(£(p',/3') + Vu(z),H r ) < d T , 

and so £(p',/3') + Vu;(z) G Sb that is, (p' + Dp(z),/3' + ^ t (z)) G db Thus (b) holds. 
The proof is now complete. □ 


5. Construction of admissible class U 

In this section, we define a suitable admissible class U as required in Section 3. 
Assume and uq are as given in Theorem 11.11 with (11. 4ft fulfilled in addition. Let 
M = II-DuoIIl^q) > 0. 
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5.1. The modified uniformly parabolic problem. We first apply Lemma [2.21 to 
construct the function / G C 3 ([0, oo)) with a suitable choice of 5 G (0,1/2) and 
1 < A < m + (5) according to the value of M as follows (see Figure 1). 

(i) If 0 < M <1, we select 0 < 6 < and arbitrary 1 < A < m + (S). 

(ii) If M > 1 and A > 0, we select 5 = i + ^J^x) 2 an< ^ arbitrary A G (M, M + A); 
in this case, m + (S) — M + A. 

Note that in both cases we have M < A. Once / G C 3 ([0, oo)) is constructed, we 
define 

A{p) = f(\p\ 2 )p (p e R n ). 

By Lemma 12.2L equation u t = di y A(Du) is uniformly parabolic. We have the follow¬ 
ing result. 

Lemma 5.1. With 5 selected above and S$ defined by (14.121) . one has 

(p, A{p)) G S$ V m_(<5) < \p\ < M. 

Proof. From the definition of set S$, it follows that (p, A(p)) G Ss if and only if 
j^| < fdP\ 2 ) < 1+ X | p | 2 ; namely, p(m_(S)) = 5 < p*{\p\) < p(\p\). 

By Lemma [2.21 this condition is satisfied if m_(<5) < \p\ < A. □ 

5.2. The suitable boundary function <F. By Theorem 12.11 the initial-Neumann 
boundary value problem 

{ vf = di v(A(Du*)) in Qt, 
du*/d n = 0 on dQ x [0, T], 

u* (x, 0) = uo(x), x E fl 

admits a unique classical solution u* G C 2+a, ^ L (Ot) satisfying 

\Du*(x,t)\ < M , (x,t) G Ot- 

From conditions HI.2)1 and fjl.4|) . we can find a function h G C 2+Q (f2) satisfying 

A h = uq in fl, dh/d n = 0 on <9fl. 

Now let v 0 = Dh G C' 1+a (0; R n ) and define, for (x,t) G f l T , 

(5.2) v*(x, t) — v 0 (x) + f A(Du*(x, s)) ds. 

Jo 

Define = (u*,v*) G C' 1 (0 T '; R 1+n ). Then it is easy to see that satisfies condition 
(11.611 above; i.e., 

{ u*(x, 0) = u 0 (x ) (x G fl), 
di vv* — u* in fl T , 

■ n| aQ = 0 V t G [0, T], 

Lemma 5.2. Let 

1C 8 = {(p,cr(p)) | \p\ < m-(S)}. 

Then 

(. Du*(x , t),v*(x, t)) G S 8 U ICg V (x, t ) G f It- 
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Proof. Given (x,t) G Cl T , let p = Du*(x,t ); then \p\ < M. By (I5.2]h vf(x,t) 
If |p| < m_(S), then A(p) = a(p) and hence 

(Du*(x,t),Vt(x,t)) = (p,A(p )) = (p,a(p)) G K s . 

If m_(5) < |p| < M, then by Lemma [5. II 

(Du*(x,t),v*(x,t)) = ( p,A(p )) G S s . 

Hence (Du*,v*) G S$ U K-s in D T . 


A(p). 


□ 


5.3. The admissible class U. In what follows, we say that a function u is piecewise 
C 1 in Q T and write u G Cp C (D T ) provided that there exists a sequence of disjoint 
open sets {Gj}JL 1 in f It with \dGj\ = 0 such that 

u G C\Gj) V j G N = {1,2, • • •}, \Q T \ U f =1 Gj\ = 0. 

(For our purpose it is also acceptable to allow only finitely many pieces in this defi¬ 
nition.) 


Definition 5.1. Let /i = ||u*||L°°(n T ) + 1- We define the admissible class 

«={»€C‘nfJ”(lV) | IklU- <k 3i>eC‘ nw^fSlriR") 

(5.3) 1 


such that div v = u and (Du, v t ) G S$ U K-s a.e. in fl T 
For each e > 0, let U e be defined by 

L/ e = juGL/|3uG Cp C D I / F,;* 0 °(Gt; R n ) such that divi; = u and 

( Du,v t ) G U JC§ a.e. in f 1 T , and J n \v t — cr(Du)\dxdt < e|fir| \- 


Remark 5.2. Clearly u* G U\ so U is nonempty. Also U is a bounded subset of 
W^^t) as Ss U Ks is bounded. Moreover, by Corollary 14.41 for each u G U, 
its corresponding vector function v satisfies ||ut||i,<x>(n T ) < 1/2. Finally, note that 
m_(<5) < \Du*\ < M on some nonempty open subset of Q.t and A(Du*) a (Du*) 
on this set; hence u* itself is not a weak solution to (11.11) . 


6 . Density property and Proof of Theorem 11.11 

In this final section, we prove the density property of the sets U e and then complete 
the proof of Theorem 11.11 

6.1. The density property of U. Let U and U e be as defined in Section 5. We 
establish the density property (13.31) . 

Theorem 6.1. For each e > 0 ,U e is dense in U under the L°°-norm. 

Proof. Let u ElA, p > 0. The goal is to show that there exists a function u ElA e such 
that \\u — w||z,°°(f 2 T ) < rj. For clarity, we divide the proof into several steps. 

1. Note || rq || £oo(q t ) < p — r 0 for some r 0 > 0 and there exists a function v G 
C; c n^°°(0 T ;R" ) such that divu = u and (Du, v t ) G Ss U/Q a.e. in fl T . Since both 
u and v are piecewise C 1 in Qt, there exists a sequence of disjoint open sets {Gj}fi =1 
in Qt with \8Gfi — 0 V) > 1 such that 

u G C 1 (G j ), v G C 1 (G j ; R n ) Vj > 1, \D T \ U f =l Gfi = 0. 
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2. Let j 6 N be fixed. Note that (Du(z), v t (z)) G S$ U K,$ for all z G Gj and 
that Hj = {z G Gj | (Du(z), v t (z)) G dS$} is a (relatively) closed set in Gj with 
measure zero. So Gj = Gj \ Hj is an open subset of Gj with | G 3 \ = \Gj\, and 
(Du(z), Vt(z)) G Ss U K$ for all z G Gj. 

3. For each r > 0, let Q r = {(p, (3) G Sg \ \/3 — cr(p)| > r, dist ((p, /3),dS$) > r)}; 
then Q t CC Ss- We can find a Tj > 0 such that 

(6.1) f \v t (z)-a(Du(z))\dz<^j\n T \, 

where z = (x,t) and Fj = {z G Gj \ (Du(z) , v t (z)) ^ Gtj}- To check this, note 
Fj = F*. U F?., where F ^ = {z G Gj \ \v t (z) — a(Du(z))\ < r} and F? = {z G 
Gj | dist((Hn(^), u f (^)), 95^) < r}. Clearly, J pl \v t — cr(Du)\dz < r|STr| —» 0 as 
r —> 0 + . Since F? is decreasing as r —* 0 + and v t = a (Du ) on fl r >oC3, it follows that 
If? \v t — a(Du)\dz —* 0 as r —* 0 + . 

4. Let Oj = Gj \ Fj. Furthermore, we may require that the number Tj be chosen 
in such a way that either Oj is empty or Oj is a nonempty open set with \dOj\ = 0 
(see [32))- Let J be the set of all j G N with Oj ^ 0. Then for j ^ J, Fj = Gj. 

5. We now fix a j G J. Note that 0 3 = {z G Gj \ (Du(z) , v t (z)) G Q Tj } and 
that Kj := Q Tj is a compact subset of Ss- Let O C R n be an open box with Q C O 

and I = (0, T). Applying Theorem 14.61 to box Q x /, Kj, and e' = we obtain a 
constant pj > 0 that satisfies the conclusion of the theorem. By the uniform continuity 
of a on compact subsets of R", we can find a 9 — 9 e j > 0 such that 

( 6 - 2 ) \ a (p) ~ a (p')\ < ^ 

if \p\, \p'\ < m + (5 ) and |p — p'\ < 6. Also by the uniform continuity of u, v, and their 
gradients on Gj, there exists a u 3 > 0 such that 

( . | u(z) - u(z ')| + |Vu(*) - Vu(z ')| + | v(z) - v(z ')| 

^ ' ’ +|Vn(^) - Vv(z')\ < min{^-, 9} 

whenever z,z' G Gj and \z — z'\ < Vj. We now cover 0 3 (up to measure zero) by a 
sequence of disjoint open cubes {Q* x F-\™d x in fl T whose sides are parallel to the 
axes with center z 3 and diameter /* < u 3 . 

6. Fix i G N and write w = (u,v), £ = ^ = Vw(Zj) = 

By the choice of pj > 0 in Step 5 via Theorem 14.61 since Q l j x I 1 - C Q x I and 
(p, fd) G Kj, for all sufficiently small p > 0, there exists a function u 1 - = (p*, A)) G 
C™(Q) x /j; R 1+n ) satisfying 

(a) div -0} = 0 in Q* x i], 

(b) (p 7 + Dp l j(z),/3' + (iifitiz)) G for all z e Q) x Pj 
and all |(p 7 ,/3 7 ) - (p,/3)| < pj, 

( c ) 11^} lll/°°(Qj x/j) < Pi 

( d ) fqixij \P + (V’i)tW - o-(p + D^j-(^))l^ < x J jl/IQ X 4 

(e) Jgi <p*(x, t)dx = 0 for all t G /*, 


(Du(Zj) u t (zfi\ 

yDv(Zj) v t (zj) J 
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(f) H0iP})t||.L° o (QjxJj) < P- 

Here, we let 0 < p < min{To, 77}, where C n is the constant in Theorem 

12.31 and C is the product of C n and the sum of the lengths of all sides of Q. By 
(e), we can apply Theorem 12.31 to pj on Qj x 7® to obtain a function r/® = TZp^ G 

x /■; R®®) fl H'o’°°(Q® x /■; R") such that div 9* = p 1 , in Q® x /■ and 


(6.4) 


\M) 


'j)t\\L°°{Q i - xPA < ^ L£ - 


Pi 

*3—3' — 2 

7. As |iy|, |ct(.Dm)| <1/2 a.e. in f It, we can select a hnite index set I C J x N so 
that 


(6.5) 

We finally define 


/ |u t (^) - cr(Z>u( 2 ))|Gh: <-|Hr|- 

^U(j.i)e(JxN)\I 

(«, 6) = K u) + XQixI^P)^) + S'}) in 
U,i)ez 


As a side remark, note here that only finitely many functions (pfi ipj+g*) are disjointly 
patched to the original (u, v ) to obtain a new function (u, v) towards the goal of the 
proof. The reason for using only finitely many pieces of gluing is due to the lack of 
control over the spatial gradients D ("0® + g*), and overcoming this difficulty is at the 
heart of this paper. 

8. Let us now check that u together with v indeed gives the desired result. By 
construction, it is clear that u G C pc fl H / /l°°(H T ), v G C pc fl W^°°(fir; R n ). By the 
choice of p in (f) as p < To, we have ||fh||L°°(n T ) < P- Next, let (j, i) G X, and observe 
that for z eQ)x J®, with (p, ft) = (Du(Zj), v t (Zj)) G G Tj , since \z — Zj\ < /* < uj, it 
follows from (16.3jl and (16.41) that 

| (Du(z),v t (z) + ( g)) t (z )) - (p,P )| < pj, 

and so (Du(z),Vt(z)) G S$, by (b) above. From (a) and div 9® = pj, for z G Q l j x J®, 
div v(z) = div(u + 0} + g})(z) = u(z) + 0 + p){z) = u(z). 

Therefore, u EU. Next, observe 



a(Du)\dz 


'UjenFj 


| v t — a(Du)\dz 


+ / \v t — a(Du)\dz + / \v t — a(Du)\dz 

^ u 0.i)e(^xN)\xQj' ^ u 0,i)ei < 3j xI j 

= h + h + 

From (16.11) and (16.51) . we have I\ + J 2 < Note also that for (j, i) G X and 

z G Q X j x J®, from (jfT3jh (j6~4|) . and (f), 

1^(2) - o-(X>it(2))| = | v t (z) + (^)t(^) + (g})t(z) - a(Du(z) + XV)(X))| 

< | v t (z) - u t (zj)| + \v t (z}) + (0])t(X) - a(Du(z}) + XV}(z))| 

+I(^)*WI + \ a ( Du ( z )) + D( P)(z)) ~ cr(Du{z) + Dp'Az))\ 
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< | + \v t (z *■) + (^)t(-) - o-(-Dw(«j) + #¥>}(*)) I 

+|ct(L>m(z]) + Dtp^z)) - cr(Du(z ) + 

Here, as (Du(z), v t {z)) G «Si, we have |£>m(z) +.Dv 3 j(^)l = — m +(^)j an d by 

flop , \Du(Zj) — Du(z )| < 6 *. From (16.21) we now have 

\a(Du(z}) + Dp){z)) - a{Du(z ) + 

Integrating the above inequality over Qb x /*, we thus obtain from (d) that 

< ||Q} x /j| + ^ X ^ 

4 12 |Q X I\ 

<\\Q)xi% 

which yields that 1% < Hence I\ +12 +13 < elfi-rl, and sow 6 W £ . Finally, from 

(c) with p < 7 ] and the definition of u, we have 

||n n||ioo(Q T ) < r). 

This completes the proof. □ 

6.2. Completion of Proof of Theorem ll.il The existence of infinitely many weak 
solutions to problem (II. ip readily follows from a combination of Remark 15.21 Theorem 
16.11 and Theorem 13.31 

To prove the last statement of Theorem 11.11 assume M = ||-Duo||l°°(Q) > 1 and 
A > 0. Recall that in this case 5 = 7 - 7 ^ 7 ^- 777 ■ By the definition of U and Corollary 14.41 
we have \Du\ < m + (S) = M + A a.e. in f \ for all u EU. On the other hand, following 
Step 3 in the proof of Theorem 13.31 every weak solution u G Q is the L °°-limit of some 
sequence Uj G U such that Duj -E Du a.e. in fl T . So these weak solutions u G Q 
must satisfy ||-Dw||l°°(q t ) < M + A. 
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